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Abstract 



The perturbative series used to extract as (Mr) from the r hadronic width exhibits slow conver- 
Cn I gence. Asymptotic Pade-approximant and Padc summation techniques provide an estimate of these 

^^ ■ unknown higher-order effects, leading to values for as (Mr) that are about 10% smaller than current 

estimates. Such a reduction improves the agreement of as (M^-) with the QCD evolution of the strong 



Q>^ ■ coupling constant from a^ (M^). 



The Particle Data Group (PDG) quotes the following values for the strong coupling constant as deter- 
mined from Z^ and r decays |l|. 

as (Mr) = 0.35 ± 0.03 (1) 

a^(Mz) = 0.119 ±0.002 (2) 

/\ • Since these determinations of ag occur at such widely separated energies, the compatibility of these 
3 . values of a^ with the QCD evolution of the coupling constant is an important test of both QCD and the 
phenomenological results used to extract the coupling constant from the experimental data. In particular, 
as (Mr) is sufficiently large that presently unknown terms from higher order perturbation theory could 
substantially alter the value of a^ (M^) extracted from the experimental data. Fade approximant methods 
provide estimates of the aggregate effect of (presently unknown) terms from higher-order perturbation 
theory ||2|, ^, ^, ^]. As shown below, the use of Fade summation to estimate such terms leads to a decrease 
in the value of ag (Afr) extracted from r decays, improving the compatibility of a^ (Mr) and ag (Mz) 
with the QCD evolution of the coupling constant. 

The QCD evolution of the coupling constant is governed by the /? function which is now known to 
4- loop order Q. Using the conventions of Q, a = — satisfies the differential equation 

T OO 

^2^=/3(a) = -a2^fta^ , « = - (3) 

11 2^ 109 38„ 2857 5033 , 325 2 

o _ 11 - 3^/ o _ 10^ - -nf _ 2 18 ^/ + 54 "/ , .^ 

/3o- ^ , /3i- , /32- ^^ (4) 



Ps = 114.23033 - 27.133944n/ + 1.5823791n^ - 5.8566958 x lO^^n) 



(5) 



Using the value as (Mz) as an initial condition, the coupling constant can be evolved to the desired energy 
using the differential equation (y). The only subtlety in this approach is the location of flavour thresholds 
where the number of effective flavour degrees of freedom Uf change. In general, matching conditions must 
be imposed at these thresholds to relate QCD with Uf quarks to an effective theory with n/ — 1 light 
quarks and a decoupled heavy quark |^]. Using the matching threshold fith defined by mq[^th) = ^^thJ 
where ruq is the running quark mass, the matching condition to three- loop order is P] 



a("^-i) {^th) = a^'^^) {^^th) 



1 + 0.1528 



("/) {^,th)] + {0.9721 - 0.0847 (n/ - 1)} [a^"/) {^lth) 



(6) 



leading to a discontinuity of as across the threshold. Thus to determine the coupling constant at energies 

(5) 

between the c quark threshold and the b quark threshold, the (3 function with n/ = 5 is used to run as 
from Mz to /^th = m-bifJ-th) = i^^b using @ as an initial condition. The matching condition (g) is then 
imposed to find the value of as (nih) which is then used as an initial condition to evolve Os to lower 
energies via the nj = 4 beta function. 

If as (Mz) is used as the input value to determine the QCD prediction of as (M,-), then one might 
legitimately be concerned about the effect of (unknown) higher-order terms in the (3 function at lower 
energies where as is larger. Pade approximations have proven their utility in determining higher-order 
terms in the /3 function. For example, using as input the four-loop (3 function in A^-component massive 
0^ scalar field theory [|l^], asymptotic Pade methods described in Section II of Q are able to predict the 
five-loop term to better than 10% of the known five-loop contributions for A^ < 4 [^, ^, |l2[. When these 
same methods are applied to QCD, the following predictions for the unknown five-loop contribution to 
the (3 function are obtained [pT|]. 






4: 
5 : 



/34 = 83.7563 
(34 = 134.56 



(7) 
(8) 



Prom these predictions, f3 functions containing [2|2] Pade approximants can be constructed to estimate the 
sum of all higher-order contributions. These Pade-summations, whose Maclaurin expansions reproduce 
/?i) /32) /^s and and the asymptotic Pade-approximant estimates (0,g) of Pi, are given by: 



m 



n 



f 
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5 : 



/3(a) 



Pia) 



25x^ 
23x^ 



1 - 5.8963a - 4.0110a^ 
1 - 7.4363a + 4.3932a2 

1 - 5.9761a - 6.9861a2 
1 - 7.2369a - 0.66390a2 



(9) 



(10) 



Thus the QCD prediction of ag {Mr) depends on only two parameters: the initial condition as (M^) 
and the position of the five-fiavour threshold defined by mb{mi,) = nib. As will be discussed below, the 
uncertainty in the Particle Data Group value M for this threshold 



4.1 GeV < mhirrib) < 4.4 GeV 



11 



has a negligible effect on the QCD prediction of Os {Mt) compared with the uncertainty in q^ (M^) (g). 

The compatibility of the experimentally/phenomenologically determined values as {Mz) and as {Mr) 

with the QCD evolution of the coupling constant can now be studied. Figure || shows the effect on as{Q) 



of progressive increases in the number of perturbative terms in the (3 function, culminating with the Pade 



summation (^,10) for (5. It is evident that the curves for as{Q) appear to converge from below to that 
generated by the Pade summation of the /3 function, since the gaps between curves of successive order 
decrease. Using the input values (^,[ll|) for the QCD evolution of q^ down from the Z^ to r mass, we 
obtain the following range of values for a^ (M^) for successive orders of perturbation theory [ p^ : 

2 - loop 0.3055 < as {Mr) < 0.3383 (12) 

3 - loop 0.3096 < as (Mr) < 0.3442 (13) 

4 - loop 0.3112 < as (M^) < 0.3468 (14) 
Pade summation 0.3119 < q^ (M^) < 0.3480 (15) 

The dominant contribution to the uncertainty (1^15) originates from q^ {Mz) — the effect of the un 



certainty in the five-flavour threshold ( [Til) is inconsequential. It is also evident that Pade-improvement 
via (^) and (^) does not alter significantly the range of as (M-r) devolving from the empirical range for 
as (Mz), as given in (g). Moreover, the ranges ( [l^ and (^) overlap the lower end of the current PDG 
range (|l]) for the extraction of as {Mr) from Rr- We note, however, that only minimal overlap occurs 
between (|l|Jll) and the previous (1996) PDG range a^ {Mr) = 0.370 ±0.033 obtained from Rr ^. This 
marginal compatibility with the RG-devolution estimates of q^ {Mr) provided the original motivation for 
us to examine the effect of (estimated) higher-order perturbative corrections on the extraction of as {Mr) 
from Rr- 

This extraction occurs by comparing the measured value of Rr with the calculated value of 5^^' , the 
purely perturbative QCD corrections to the tree diagram for the calculation of Rr'- 



Rr ^ ^i-^-^ + '-^'^^^^^) = 3.058 

T {t ^ Ur + e + Ue) 



1 + (5(°) - (0.014 ± 0.005)1 . (16) 



The above expression, as quoted from the current PDG |13|, is based on the seminal analysis of Braaten 



Narison, and Pich |14], as well as more recent work by Neubert |15]. In particular, the numerical factor 
in parentheses represents non-perturbative contributions that are dominated by dimension-6 terms. In 
the MS scheme, the purely perturbative QCD contributions to Rr are 

l + ,5(°) = l + a(M^) + 5.2023[a(M^)]V26.366[a(M^)]^ . (17) 

where a {Mr) = as {Mr) /vr. Using the empirically motivated test value a^ {Mr) = 0.3500, corresponding 
to the central value of the PDG range (|l|), we find that the contributions of successive terms in (^) to 
6^^' are respectively 

(5^°) = 0.1114 + 0.06457 + 0.03646 (18) 

illustrating the slow convergence of perturbative field theory at the mass scale n = Mr- The ratio of 
successive terms within 6^^' appears to be nearly 0.6, indicative of significant further contributions from 
corrections to ( [T^ ) beyond three- loop-order. 

Asymptotic Pade-approximant methods |^ ^ can be utilized to estimate the aggregate effect of higher 
order terms in the series (p!7|). Given a field-theoretical perturbative series of the form 

S = 1 + Ria + R2a^ + Rsa^ + R^a"^ (19) 

in which the coefficients Rk are characterized by asymptotic behaviour R^ ~ k\C k^ [16|, the coefficient 



-R4, which we assume to be unknown, can be estimated from the known terms i?i, R2, and -R3. An [A^|M] 



Pade-approximant for the series ( p!9|) then yields coefficients R^ that differ from those of the series 
itself via the asymptotic error formula ^, ^ 

Rn+m+1 [N + M{l + c) + h]^ ' 

where {A,c,b} are constants to be determined. For example, a [2|1] Pade approximant to (^) would 
lead to the prediction 

^pade ^ |3 (21) 

It has been shown elsewhere [O] that the error formula ( pO| ) can be utilized in conjunction with [0|1] 
and [1[1] approximants to determine A and (c + 6), thereby leading to the following "asymptotic Pade- 
approximant" (APAP) estimate for R^: 

vAPAP ^ -R3 [3 + (c + h)] ^ i?3 [R2 + R1R2R3 - ^fiifis] ^22) 

^ iia [3 + (c + 6) - A] R2 [2R^ - RfR^ - RfR^] ^ ' 

If just the [2|1] Pade-approximant is used to estimate [via (pi|)] the af contribution to 5^^\ it is found 
that II 

1 + 5(0) = 1 + (M^) + 5.2023 [a (M^)]^ + 26.366 [a (M^)]^ + 109.2 [a (M^)]^ (23) 

an estimate very close to that obtained by Kataev and Starshenko using other methods [17|. The APAP- 



estimate of the a^ term, obtained via (^2|), is also positive and somewhat (20%) larger: 

1 + 5(0) = 1 + {Mr) + 5.2023 [a {Mr)f + 26.366 [a {Mr)f + 132.44 [a (M^)]^ (24) 

It is significant to note that this prediction is very close to the maximum estimated size of the fourth 



order effect used to determine the theoretical uncertainty in |14], indicating an underestimate of the 
higher order effects in the extraction of ag (-^r) from R^-. However, even if 5^^' includes an estimate of 
the [a {Mr)] term, one sees from (24) that the convergence of the perturbative series remains too slow 
to justify a truncation. For example, if Ug {Mr) = 0.3500, then the contribution of successive orders to 
5^^' is seen to be 

(5(°) = 0.1114 + 0.06457 + 0.03646 + 0.02040 . (25) 

The (estimated) fourth term is 18% of the leading perturbative contribution. Such slow convergence 
indicates that further higher order terms should contribute significantly to 5^^'. A Pade-summation, 
in this case the [2|2] approximant whose first five Maclaurin expansion terms replicate the series (24), 



provides an estimate of the total effect of higher order terms in a perturbation series ||2[. This Pade 
summation is given by 

^ ^ ^(0) ^ 1 - 6.54830 {Mr) + 10.5030 [a (MQ]' 
1 - 7.5483a {Mr) + 12.8514 [a {Mr)f 



Figure ^ compares the dependence of 6^^' on a^ {Mr) obtained from (17, 23, ^, 26). These curves 
correspond respectively to 



Truncation: Truncation of contributions to 6^^' beyond three-loop order (IT 



[2|1]: Inclusion via (^) of the [2|1] Pade-approximant estimate of the four-loop contribution to 



• APAP: Inclusion via (24) of the asymptotic error- formula estimate of the four- loop contribution 

to 5(0); 

• PS: Fade-summation estimate ( Pq) of all higher-loop contributions to 6^^'. 

Near the PDG value a (M,-) = 0.350, the Fade effects lead to a significant increase in 6^^' . In compar- 



ison with the three-loop perturbative result (17), the size of this enhancement obtained from the Fade 
summation ( p6| ) is roughly twice the enhancement obtained by including Fade estimates of the four-loop 
contributions (^,^), indicating the significance of the higher-order effects estimated in the [2|2] Fade 
summation. 

In Table ll|, we display the values for Os (Mr) one obtains for a given value of 6^^' by inverting equations 



(|17D , (|23D, (12 4|) or (|26D . The present empirical range Rr = 3.642 it 0.024 y, can be used in conjunction 



with (^) to extract the following range for the purely-perturbative correction 6^^' : 

(5^°) = 0.2048 ± 0.0129 (27) 

Using Table |l|, we find that this empirical range for 6^^' determines a corresponding range for a^ (M-r) 
for each case listed above :|j 

Truncation : Os (Mr) = 0.342 ± 0.013 (28) 

[2|1] : as {Mr) = 0.329 ±0.011 (29) 

AFAF : as (Mr) = 0.326 ± 0.011 (30) 

FS : as (Mr) = 0.314 ± 0.010 (31) 

It is evident from the above results that progressively sophisticated Fade-estimates of higher-order cor- 
rections to 5^^^ lead to progressively lower values for ag (Mt-).F| In particular, the ranges (p9|) and (pO) 
are fully enclosed within the ranges (|l^ ) and (15) from RG-devolution from as (Mz)- The range (|3l|). 



which we regard as the most accurate reflection of cumulative higher order corrections, is almost entirely 
enclosed by the lower portion of these RG-ranges. By contrast, only the lower half of the FDG "trunca- 
tion range" (^) for Os (Mt) is in agreement with the RG-ranges ( [l^ ) and (p!5|), although the range quoted 
in (p^) does somewhat better than this. 

We therefore conclude that higher-order corrections to 6^^' appear to lower the value of ag (M-r) ex- 
tracted from Rr by approximately 10%, but that this lowering seems to improve the overall compatibility 
of as {Mr) with the QCD evolution of Os from the present empirical range for as {Mz)- Alternatively, 
one can conclude that the theoretical uncertainty in as (Mr) associated with truncation of contributions 
to jC^-* past three-loop order is not likely to be bi-directional, as indicated in Q, but is rather an 0(10%) 
effect in the downward direction. 

Acknowledgment: We are grateful for financial support from the Natural Sciences and Engineering 
Research Council of Canada. 



^The central value for "truncation" increases to as (Mr) = 0.347, consistent with the central value in flj, provided we 
ilize directly the expressions given in |l4| for the non-perturbative contributions, which are weakly Qs dependent. 
^A range similar to (E9|) will also follow from the estimate of the ai correction to S^°' given in jlTJ. 
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Table 1: Values of q^ (Mr) for given values of 5^°) obtained by inverting (|T^ ["Truncation], (|2|) ["[2|1]"], 
"^ ["APAP"], and @ ["PS"] as discussed in the text. 
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Figure 1: Effect of increasing the order of perturbation theory in the QCD evolution of the strong 
couphng constant, using as (Mz) as an initial condition. Higher-loop terms in the f3 function progressively 
increase ag from the 2-loop order bottom (solid) curve to the Fade summation top (dashed-dotted) curve, 
sandwiching the three- and four- loop curves. 




Figure 2: Values for 6^^' as a function of Ug (Mr) using the four different treatments of 6^^' . The sohd 
curve uses ([l7| ) ["Truncation], the dotted curve uses (23) ["[2|1]"], the dashed-dotted curve uses (24) 
["APAP"], and the dashed curve uses (|2|) ["PS"], as discussed in the text. 



